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Abstract. The recent developments in studies of TeV radi¬ 
ation from blazars are highlighted and the implications of 
these results for derivation of cosmologically important in¬ 
formation about the cosmic infrared background radiation 
are discussed. 


1 Introduction 

Since early 90’s many papers have been published with an 
unusual combination of two keywords - Blazars and Cosmic- 
Infrared-Background (CIB) radiation. Formally, one may ar¬ 
gue that there is no apparent link between these two topics. 

• The blazars constitute a sub-class of AGN dominated by 
highly variable (several hours or less) components of broad¬ 
band (from radio to gamma-rays) non-thermal emission pro¬ 
duced in relativistic jets pointing close to the line of sight. 

• CIB is a part of the overall diffuse extragalactic background 
radiation (DEBRA) dominated by thermal emission compo¬ 
nents produced by stars and dust, and accumulated over the 
entire history of the Universe. 

While the blazars may serve as ideal laboratories for study 
of MHD structures and particle acceleration processes in rel¬ 
ativistic jets, CIB carries crucial cosmological information 
about the formation epochs and history of evolution of galax¬ 
ies. To a large extent, these two topics are relevant to quite 
independent areas of modern astrophysics and cosmology. 
Yet, the current studies of CIB and blazars, more specifically 
the sub-population of blazars emitting TeV gamma-rays (TeV 
blazars), are tightly coupled through the intergalactic (IG) 
absorption of TeV radiation by infrared photons of DEBRA. 

The astrophysical/cosmological importance of this inter¬ 
esting effect (Nikishov, 1962; Gould and Schreder, 1966; 
Jelly, 1966; Stecker et ah, 1992) was clearly recognized af¬ 
ter the discovery of TeV 7 -rays from two BL Lac objects - 
Mkn 421 and Mkn 501 (for review see e.g. Vassiliev 2000). 
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2 Cosmic Infrared Background Radiation 

CIB basically consists of two emission components produced 
by stars and partly absorbed/re-emitted by dust during the 
entire history of evolution of galaxies. Consequently, two 
distinct bumps in the spectral energy distribution (SED) of 
red-shifted radiation at near infrared (NIR) A ~ 1-2 pm and 
far infrared (FIR) A ~ 100-200 //in wavelengths, and a mid 
infrared (MIR) “valley” between these bumps are expected. 
Because of the heavy contamination caused by foregrounds 
of different origin, predominantly by the zodiacal (interplan¬ 
etary dust) light, the measurements of CIB contain large un¬ 
certainties. Moreover, these results only conditionally can 
be treated as direct measurements, because their interpreta¬ 
tion primarily depends on the modeling and removal of these 
foregrounds. Therefore, the direct observations of CIB gen¬ 
erally allow derivation of the flux upper limits rather than 
detection of positive residual signals. In Fig. 1 we show the 
CIB fluxes based on the latest reports, and refer the reader to 
the review article by Hauser and Dwek (2001) on the current 
status of direct observations of CIB. 

Presently, the most reliable results are obtained from the 
COBE observations at NIR and FIR domains where the con¬ 
tribution of the zodiacal light becomes comparable with the 
CIB flux. This concerns, first of all, the FIR wavelengths at 
140 and 240 pm (Hauser et al., 1998, Schlegel et al., 1998, 
Lagache et al., 1999), and perhaps also at 100 pm (Lagache 
et al., 1999; Finkbeiner et al., 2000) where the results of 

3 different groups are in a reasonable agreement with each 
other, but somewhat higher compared to the theoretical pre¬ 
dictions. If the reference of these fluxes to the truly extra- 
galactic background is correct, this would imply that most of 
the star formation in the early Universe must occur in highly 
obscured, dusty environments. 

Detections of CIB are claimed also at NIR - at 2.2 pm and 
3.5 /im wavelengths (Dwek and Arendt, 1998; Gorjian et 
al., 2000; Wright and Johnson, 2001). These fluxes together 
with the J-band upper limit at 1.25 pm (Wright and John¬ 
son, 2001 ) and the fluxes of the optical light derived from the 
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Fig. 1 . SED of Cosmic Background Radiation. The reported fluxes 
at 140 and 240 gm (stars) are from Hauser et al. (1998), at 100 /im 
(filled triangle) from Lagache et al. (1999), at 60 and 100 /im (open 
circles) from Finkbeiner et al. (2000). The point at 15 (im (open 
square) is derived from the ISOCAM source counts (Franceschini et 
al., 2001). The fluxes reported at NIR and optical wavelengths are 
shown by filled squares (Wright and Johnson, 2001), open circles 
(Matsumoto, 2000), open diamonds (Cambresy et al., 2001), and 
filled circles (Bernstein, 1998). The open triangles from Pozzetti et 
al. (1998) correspond to lower limits. While the reference models 
2 and 3 should be considered as lower and upper limits, the curve 1 
may be treated as the “model of choice”. 


HST data below 0.8 /mi (Bernstein, 1998), agree with the 
recent theoretical calculations by Primack et al. (2001). An 
independent analysis of the COBE data by Cambresy et al. 
(2001), as well as the results of the Japanese IRTS satellite 
(Matsumoto et al. 2001) are in agreement with the fluxes at 
2.2 /rm and 3.5 /tm reported by Wright and Johnson (2001), 
but at shorter wavelengths the claimed fluxes are noticeably 
higher compared to other measurements (see Fig. 1). 

The “best guess” estimate of the SED at optical/NIR wave¬ 
lengths of about 20-50 nW/m 2 sr is comparable with the FIR 
flux of about 40-50 nW/m 2 sr (Madau and Pozzetti, 2000). 
This indicates that an essential part of the energy radiated 
by stars is absorbed and re-emitted by cold dust. Our cur¬ 
rent knowledge of MIR, which carries information about the 
warm dust component, is quite limited. The only flux esti¬ 
mates in this band derived from the ISOCAM source counts 
at 6 and 15 /im (Franceschini et al., 2001) apparently should 
be taken as lower limits. 

These measurements stimulated new calculations of CIB 
based on different phenomenological and theoretical (back¬ 
ward, forward, cosmic chemical evolution, semi-analytical, 
etc.) approaches. In spite of certain achievements, the abil¬ 
ity of current (most successful) models to reproduce the CIB 
measurements nevertheless should not be overstated, because 
all these models contain a number of adjustable parameters 
and in fact are “primarily designed for that purpose” (Hauser 
and Dwek, 2001). 

In Fig. 1 three CIB models are shown which provisionally 



Fig. 2. Mean free-path of VHE 7 -rays in the intergalactic medium 
calculated for 3 CIB models presented in Fig. 1. 


may be assigned as “nominal” ( 1 ), “low” ( 2 ), and “extreme” 
(3) descriptions of CIB. The curve 2 corresponds to the so- 
called “LCDM-Salpeter” model of Primack et al. (1999), 
but assuming twice larger contribution of the dust compo¬ 
nent. Similar NIR spectra with a flux of about 10 nW/m 2 sr 
at 1 /rm have been assumed by De Jager and Stecker (2001) 
and Franceschini et al. (2001). The curve 1 is close to the 
prediction by Primack et al. (2001) based on the so-called 
“Kennicut” stellar initial mass function (IMF). And finally, 
the curve 3 is “designed” to match the extreme fluxes re¬ 
ported by Matsumoto (2000) and Cambresy et al. (2001) at 
NIR and Finkbeiner et al. (2000) at FIR. 


3 Absorption of TeV gamma-rays in CIB 


To calculate the mean free path of 7 -rays A (E) in the IGM 
one must convolve the CIB photon number distribution n(e) 
with the pair production cross-section. Because of the nar¬ 
rowness of the latter, for broad-band photon spectra over half 
the interactions of a 7 -ray photon of energy E occur with a 
quite narrow interval of target photons AA ~ (1 ± 1/2)A* 
centered on A* ~ 1.5(E/1 TeV) /im. This gives a conve¬ 
nient approximation for the optical depth t{E ) = d/A(E) 
for 7 -rays emitted by a source at a distance d: 


t(E) ~ 1 


mcib(A*) \ 
10 nW/m 2 sr J 



where «cib = vl(v) = e 2 n(e) is the SED of CIB, z is 
the source redshift, Hqq is the Hubble constant normalized 
to 60 km/s Mpc. No deviation of the observed 7 -ray spec¬ 
trum from the intrinsic (source) spectrum, e.g. by a factor 
of < 2, would imply t(E) < In 2, and consequently pro¬ 
vide an upper limit on the CIB flux at A *(E). For example, 
this condition for 1 TeV 7 -rays from Mkn 421 or Mkn 501 
(z — 0.03) gives ucib < 20 nW/m 2 sr around 1-2 /im. 
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Fig. 3. Absorption-corrected SED of Mkn 501 reconstructed for 
the analytical presentation of the measured spectral points given by 
equation (a) and assuming 3 different CIB models shown in Fig. 1. 
The experimental points (filled circles) correspond to the 1997 time- 
averaged spectrum of Mkn 501 (Aharonian et al., 1999a). The star 
corresponds to the recently re-estimated flux at E « 21 TeV (Aha¬ 
ronian et al., 2001a). The vertical line at 17 TeV indicates the edge 
of the spectrum measured with high statistical significance. 


Since this estimate is quite close to the recent NIR measure¬ 
ments (see Fig. 1), we may conclude that already at 1 TeV we 
detect significantly absorbed 7 -radiation from Mkn 421 and 
Mkn 501. This is demonstrated in Fig. 2 where we present 
accurate numerical calculations of the mean free path of 7 - 
rays for three CIB reference models shown in Fig. 1. The 
horizontal lines indicate the distances to the reported TeV 
blazars (see below), as well as to 2 nearby prominent extra- 
galactic objects - the radiogalaxy M 87 and the ultralumi- 
nous starburst galaxy Arp 220. The curve 2 calculated for 
the CIB lower limit (model 2 in Fig. 2), should be treated as 
an upper limit for the mean free path. This implies that we 
cannot ignore the IG-absorption of TeV 7 -rays from sources 
located beyond 100 Mpc. On the other hand, this gives a 
unique chance to extract information about CIB by detecting 
absorption features in TeV spectra of extragalactic objects. 

Very often such a possibility is reduced to the search for 
sharp cutoffs in the energy spectra of extragalactic TeV sources 
(see e.g. Stecker et al., 1992). However, in many cases this 
could be a misleading recommendation. Indeed, the cutoff 
in a 7 -ray spectrum does not yet imply an evidence for the 
IG-absorption and vice versa, the lack of the cutoff cannot 
be interpreted as absence of IG-absorption. It is interest¬ 
ing to note that all CIB models with conventional spectral 
shape between 1 and 10 pm, predict almost constant (energy- 
independent) mean free path of 7 -rays, and correspondingly 
insignificant spectral deformation at energies between 1 and 
several TeV (Primack et al., 2001). The explanation of this 
effect is straightforward. If the spectrum of background pho¬ 
tons in a certain energy interval has a power-law dependence. 


n(e) oc e -3 , the mean free path in the corresponding 7 -ray 
energy interval A (E) oc i? 1- - 3 . Within \pm < A < 10 p 
the CIB spectrum is approximately described by a power- 
law n(e) oc e -1 or wcib oc A -1 (see Fig. 1), therefore in the 
interval between 1 and several TeV the 7 -ray mean free path 
only slightly depends on energy (Fig. 2). In order to demon¬ 
strate the impact of the IG-absorption, in Fig. 3 we show 
the SED of Mkn 501 together with the absorption-corrected 
spectra reconstructed for 3 CIB models shown in Fig. 1. 

4 Observations of TeV Blazars 

Many nonthermal extragalactic objects representing differ¬ 
ent classes of AGN and located within 1 Gpc are potential 
TeV sources. First of all this concerns the BL Lac popula¬ 
tion of blazars of which two nearby representatives, Mkn 421 
(z = 0.031) and Mkn 501 (z = 0.034) are firmly established 
as TeV 7 -ray emitters. The current list of extragalactic TeV 
sources contains 4 more BL Lac objects with reported signals 
at 4 to 6 sigma level: 1ES 2344+51 [Whipple (Catanese et al., 
1998)] at z = 0.044, PKS 2155-304 [Durham (Chadwick et 
al., 1999)] at z = 0.116, 1ES 1959+650 [Telescope Array 
(Nishiyama et al., 1999) and HEGRA (Gotting et al., 2001)] 
at 2 = 0.048, and 1ES 1426+428 [VERITAS (Horan et al., 
2001) and HEGRA (Gotting et al., 2001)] at z = 0.129. 

Since the discovery of TeV 7 -radiation of Mkn 421 (Punch 
et al. 1992), this object has been subject of intensive studies 
through multiwavelength observations. The TeV flux of the 
source is variable with typical average value between 30% 
to 50% that of the Crab Nebula, but with strong and rapid, 
as short as 0.5 h, flares (Gaidos et al., 1996) which corre¬ 
late with the source activity at other wavelengths. Until the 
spectacular high state of the source in 2001 , the spectral stud¬ 
ies were based mainly on the data taken during a quiescent 
or moderately high states. In particular, the spectrum mea¬ 
sured by HEGRA during the 1998 “ASCA” campaign (Taka- 
hashi et al., 2000 ) is fitted by a steep power-law with photon 
index T ~ 3 (Aharonian et al., 1999b). The exceptionally 
bright and long-lasting activity of Mkn 421 in 2001 allowed 
the VERITAS (Krennrich et al., 2001) and HEGRA groups 
(Horns et al., 2001) to derive the time-averaged gamma-ray 
spectrum of the source in the high state which up to ~ 15 
TeV is described as dN/dE = I\E r exp (—E/E 0 ), i.e. 
by the same canonical “power-law with exponential cutoff” 
function found earlier for the high-state, time-averaged spec¬ 
trum of Mkn 501 (Aharonian et al., 1999a). The spectra 
of these sources are not, however, identical. In Fig. 4 we 
show the energy spectra based on approximately 40,000 and 
60,000 TeV 7 -ray events detected by HEGRA during the ex¬ 
ceptionally high states of Mkn 501 in 1997 and Mkn 421 in 
2001, respectively. They are described by different [F, Eq] 
combinations: T = 1.92 ±0.03 and Eq = 6.2 ±0.4 TeV for 
Mkn 501 (Aharonian et al., 1999a) and T = 2.23 ± 0.04 and 
Eq = 4.0. ± 0.4 TeV (Horns et al., 2001). The power-law 
part of the spectrum of Mkn 421 is steeper, and the exponen¬ 
tial cutoff starts earlier. Since both sources are located at ap- 
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Fig. 4. Energy spectra of Mkn 501 (Aharonian et al., 1999a) and 
Mkn 421 (Horns et al., 2001) as measured by HEGRA CT system 
during the high states of sources in 1997 and 2001, respectively. 


proximately same distances, the difference in the cutoff ener¬ 
gies may be interpreted as an evidence against the hypothesis 
which attributes the cutoffs to the pure IG-absorption effect. 

For the typical energy resolution of Cherenkov telescopes 
of about 20 % or larger, one has to be careful with conclu¬ 
sions concerning the spectral shape at energies beyond « 

3 /%. Obviously for determination of spectra with sharper 
(super-exponential) cutoffs, a significantly better energy res¬ 
olution is required. Motivated by this, HEGRA collabora¬ 
tion recently re-analyzed the “1997 high-state” spectrum of 
Mkn 501 using an improved method for energy reconstruc¬ 
tion with ~ 12 % resolution achievable for the stereoscopic 
mode of observation (Hofmann, 2000). The new analysis 
confirms the result of the previous study up to 17 TeV, but 
constrains stronger the flux at 21 TeV (Aharonian et al., 2001a) 
With this new point, the Mkn 501 spectrum at very high en¬ 
ergies is better described by a super-exponential cutoff (com¬ 
pare curves (a) and (b) in Fig. 3). 

Despite the huge overall photon statistics, the spectra of 
both sources above 10 TeV can be determined using data ac¬ 
cumulated over long periods of observations. Obviously, the 
time-averaged spectra of variable sources obtained in such 
way can be considered as astrophysic ally meaningful pro¬ 
vided that the spectral shape is essentially time-independent, 
i.e. does not correlate with the absolute flux. Remarkably, 
it turned out that during the high state in 1997 the shapes of 
daily spectra of Mkn 501remained essentially stable despite 
dramatic flux variations (Aharonian et al., 1999a), although 
for some specific flares significant spectral variations cannot 
be excluded. For example, the CAT group (Djannati-Atai et 
al., 1999) found an evidence for a non-negligible hardness- 
intensity correlation for the 1997 April 7 and 16 flares. 

The sensitivity of current ground-based detectors is not 


o 
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Time (sec from T 0 = 11987.0 [TJD]) 

Fig. 5. Simultaneous X-ray and 7 -ray observations of the 2001 
March 19 flare by RXTE (squares) and VERITAS (circles). The 
bottom panel shows the RXTE/PCA “6-11 keV/3.5-4.5 keV” hard¬ 
ness ratio (from Fossati et al., 2001). 


sufficient to study the spectral variability of TeV radiation of 
Mkn 421 and Mkn 501 in quiescent states on timescales less 
than several days. However, even the time-averaged spec¬ 
tra of these sources in low states are of certain interest. The 
HEGRA observations of Mkn 501 in 1998 and 1999, when 
the source was at a ~ 10 times lower flux level compared 
to the average flux in 1997, revealed a noticeable steepening 
(by 0.44 ± 0.1 in photon index) of the energy spectrum (Aha¬ 
ronian et al. 2001b). This seems to be true also for Mkn 421, 
al least at low energies. While the spectrum of this source in 
a quiescent state has a steep power-law spectrum with pho¬ 
ton index T ~ 3, the time averaged spectrum of the 2001 
outburst is significantly flatter, with T ~ 2.2 around 1 TeV. 

The SEDs of BL Lacs are expected to be very hard up to 
TeV energies, therefore the imaging atmospheric Cherenkov 
telescopes (IACTs) are nicely suited to search for short sig¬ 
nals from these objects. In particular, the HEGRA, CAT 
and VERITAS IACTs can follow m 10 -11 erg/cm 2 s (at 1 
TeV) flares of Mkn 421 and Mkn 501 on timescales less than 
several hours, and thus are well-matched to the sensitivity 
and spectral coverage of X-ray satellites like RXTE, Bep- 
poSAX and XMM for multiwavelength monitoring of flux 
variations. This is a key condition which makes the simulta¬ 
neous X- and 7 -ray observations meaningful and very impor¬ 
tant (e.g. Pian et al., 1998, Sambruna et al., 2000; Maraschi 
et al., 1999; Kataoka et al., 1999; Rrawczynski et al., 2000; 
Takahashi et al., 2000). The best so far results in this regard 
became available recently, after the well coordinated multi¬ 
wavelength campaigns of spectacular flares of Mkn 421 in 
2001. On several occasions, truly simultaneous observations 
by RXTE and TeV instruments with durations up to 6 hours 
per night were carried out. A nice sample of such events. 
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Fig. 6. Time-evolution of the absolute flux and the hardness ratio of 
TeV radiation of Mkn 421 observed by the HEGRA I ACT system 
on March 21/22, 2001 (Homs et al., 2001). 


the 2001 March 22 flare detected by RXTE and VERITAS 
(Fossati et al., 2001) shows a clear keV/TeV correlation on 
sub-hour timescales (see Fig. 5). 

For deep understanding of acceleration and radiation pro¬ 
cesses in jets, it is crucial to search for spectral variability on 
timescales comparable to the characteristic dynamical times 
of about 1 h. The HEGRA IACT system is able to perform 
such studies, if the flux at 1 TeV exceeds 1-2 Crab units, 
Jcn&b{E > 1 TeV) ~ 1.7 x 10~ n ph/cm 2 s. The strongest 
flares of Mkn 421 in 2001 provide us with such unique data. 
The preliminary spectral analysis of observations of these 
flares (Horns et al., 2001) shows a complex picture. If the 
spectra of some of these flares do not show noticeable spec¬ 
tral changes, the others, in particular the 21/22 March 2001 
flare demonstrates an impressive correlation of the hardness 
ratio with the absolute flux (see Fig. 6 ). 

5 Radiation mechanisms of TeV gamma-rays 

The studies of spectral evolution of TeV 7 -rays and their cor¬ 
relations with X-rays on sub-hour timescales not only im¬ 
ply a new standard in 7 -ray astronomy, but also indicate that 
the quality of the TeV data is approaching to the level which 
would allow us to follow and resolve simultaneously the pre¬ 
dicted fluctuations in the putative synchrotron and IC emis¬ 
sion components on timescales close to the shortest ones likely 
in these objects. This may have two important consequences 
(Coppi and Aharonian, 1999a): 

• Matching the observed X-ray/TeV light curves (as opposed 
to simply fitting snapshot spectra obtained many hours and 
days apart) should provide a very stringent test of the so- 



Fig. 7. Interpretation of the quasi-simultaneous X-ray (BeppoSAX) 
and TeV (CAT) observation of 1997 April 7 and 16 flares of Mkn 
501 within the homogeneous SSC model (Tavecchio et al., 2001). 


called synchrotron-self Compton (SSC) model since we have 
two detailed handles on the single electron distribution re¬ 
sponsible for both emission components. This test can rule 
out alternative hadronic models which are less attractive but 
still viable options for explanation of TeV emission. 

• If the SSC model works, it will be possible to fix the key 
model parameters and calculate the blazar’s intrinsic spec¬ 
trum, comparing the observed variations of absolute duxes 
and spectral shapes with the predictions of self-consistent, 
time-dependent numerical codes. This is a crucial point be¬ 
cause with an estimate of the intrinsic TeV spectrum, then, 
and only then, we can attempt to estimate the IG-absorption 
effect by comparing the intrinsic and observed spectra, and 
thus to get information about CIB. 

The observed X/TeV correlations are often interpreted as 
a strong argument in favor of SSC model, the simplified one- 
zone version of which generally gives satisfactory fits to the 
observed X-ray and TeV spectra of Mkn 421 and Mkn 501 
(e.g. Mastichiadis and Kirk, 1997; Bednarek and Protheroe, 
1999; Kataoka et al., 1999; Krawczynski et al., 2000; Tavec¬ 
chio et al., 2001, etc.). An example of very good spectral 
fits for the 1997 April flares of Mkn 501, obtained for fixed 
“standard” model parameters <5j = 10 and B ~ 0.3 G and 
changing only the maximum energy of accelerated electrons, 
is shown in Fig. 7 (Tavecchio et al., 2001). It should be no¬ 
ticed, however, that the calculation in Fig. 7 do not take into 
account the IG-absorption of 7 -rays which otherwise results 
in rather flat intrinsic 7 -ray spectra (Coppi and Aharonian, 
1999b; Konopelko et al., 1999; Guy et al, 2000). In particu¬ 
lar, any CIB model that predicts NIR fluxes close or higher 
the CIB reference model no.l would require unusually hard 
intrinsic TeV spectrum (Guy et al., 2000) with photon in¬ 
dex T < 1.5 (see Fig. 3). Such spectra cannot be easily 
explained by the standard one-zone SSC model without in¬ 
voking extreme jet parameters like Doppler factors <5j ~ 100, 
and very small B-field B < 0.01 G, which however reduce 









6 


the radiation efficiency to an almost unacceptably low level 
(Krawczynski et al., 2001). The situation is even more dif¬ 
ficult for Mkn 421, the X-ray synchrotron cutoff of which 
never extends beyond 10-20 keV. This significantly limits the 
freedom for assuming parameters which would allow hard 
multi-TeV IC spectra. 

Despite these difficulties, the leptonic models (perhaps in 
their more sophisticated versions including the multi-zone 
SSC and/or External Compton components) remain the gen¬ 
erally preferred concept for TeV blazars. These models have 
two attractive features: (i) the capability of the (relatively) 
well understood shocks to accelerate electrons to TeV ener¬ 
gies (Sikora and Madejski, 2001; Pelletier, 2001) and (ii) ef¬ 
fective production of tightly correlated X-ray and TeV 7 -ray 
emission components via the synchrotron and inverse Comp¬ 
ton radiation channels (e.g. Ulrich et al., 1997). 

The so-called hadronic models are generally lacking in 
these virtues. These models assume that the observed 7 -ray 
emission is initiated by accelerated protons interacting with 
ambient mater (Bednarek, 1993; Dar and Laor, 1997; Pohl 
and Schlickeser, 2000), photon fields (Mannheim, 1996), mag¬ 
netic fields (Aharonian, 2000) or both (Miicke and Protheroe, 
2001). The models involving interactions of protons with 
photon and B-fields require particle acceleration to extreme 
energies exceeding 10 19 eV which is possible if the accelera¬ 
tion time is close to f acc = t?(r g /c) with 77 ~ 1 (r g and 77 are 
the so-called gyro-radius, and gyro-factor, respectively, c is 
the speed of light). This corresponds (independent of specific 
acceleration mechanism) to the maximum (theoretically pos¬ 
sible) acceleration rates which hardly can be achieved by the 
conventional diffusive shock acceleration mechanism (e.g. 
Pelletier, 2001, Sikora and Madejski, 2001). On the other 
hand, the undisputed (for any TeV blazar model) condition 
of high (close to 100 %) efficiency of radiative cooling of ac¬ 
celerated particles requires extreme parameters characteriz¬ 
ing the sub-parsec jets and their environments, in particular 
very high (often, almost unrealistic) densities of the thermal 
plasma, radiation and/or B-fields. But of course, these prob¬ 
lems should not prevent us from considering the hadronic 
models as a viable alternative to the leptonic models. 

It should be noticed in this regard that the very fact of 
strong X/TeV correlations do not yet exclude the hadronic 
models, because the secondary electrons - the products of 
interactions of accelerated protons and primary 7 -rays - can 
be readily responsible for the synchrotron X-radiation. Also, 
some of the hadronic models, for example the synchrotron- 
proton model can provide not only effective radiative (syn¬ 
chrotron) cooling but also quite good TeV spectral fits. The 
proton-synchrotron model allows very hard intrinsic TeV spec¬ 
tra and can rather naturally explain the stable spectral shape 
of TeV radiation of Mkn 501 during strong flares by the self- 
regulated synchrotron cutoff at i5 cut ~ 0 . 3 r 7 — 1 (5j (Aharonian, 
2000). Fig. 8 demonstrates the ability of the synchrotron- 
proton model to explain the spectra of TeV radiation of Mkn 
501 in both the low and high states of the source. 

In the proton-synchrotron model we deal with highly mag¬ 
netized (100 G or so) condensations of 7 -ray emitting clouds 



Fig. 8. Interpretation of TeV radiation of Mkn 501 in high and 
low states by the proton-synchrotron model. The dashed and solid 
curves correspond to the spectra before and after correction for IG- 
absorption, respectively (Aharonian, 2000). 


of EHE protons, where the magnetic pressure dominates over 
the pressure of relativistic protons. In the SSC models the sit¬ 
uation is exactly opposite. The pressure of relativistic elec¬ 
trons is much larger than the pressure of the B-field (Aharo¬ 
nian, 2000; Kino et al., 2001). Surprisingly, for TeV blazars 
we do not have an intermediate arrangement! 

Both the leptonic and hadronic models of TeV blazars pose 
a number of questions. The answer to some of these question 
may become possible in the near future. Such an optimistic 
view is based on the recent progress in the theoretical (time- 
dependent) studies of characteristics of nonthermal radiation 
produced in AGN jets, and, of course, on the unique data ob¬ 
tained by RXTE and ground-based 7 -ray detectors during the 
spectacular flares of Mkn 421 in 2001. Importantly, on some 
occasions positive signals from Mkn 421 were detected also 
by low-energy air Cherenkov instruments STACEE (Ong et 
al., 2001) and CELESTE (Le Gallou et al., 2001). These re¬ 
sults are obtained at energies around 100 GeV at which the 
IG-absorption could be safely neglected, and therefore the re¬ 
ported fluxes can be used as “calibration” points when com¬ 
paring the theoretical predictions with experimental data. 

Because Mkn 421 and Mkn 501 are located at almost same 
distances, we can require that any spectral feature attributed 
to the IG-absorption be exactly the same for both sources. 
This is an important circumstance which can help to dis¬ 
entangle the local features in the intrinsic spectra from the 
contribution induced by interactions with the CIB photons. 
At the same time it is crucial to have TeV blazars at larger 
distances, so we can exploit the differences in the spectral 
modification factors. Fortunately, in addition to Mkn 421 
and Mkn 501, presently we do have 4 more TeV blazar can¬ 
didates - 1ES 2344+5, 1ES 1959+65, PKS 2155-304 and 
1ES 1426+428 - located at larger distances (see Fig. 2). 
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Fig. 9. Impact of the intergalactic absorption on the TeV spectra 
of Mkn 501 (upper curves) and 1ES 1426+428 (lower curves). The 
curves marked as “1”, “2”, and “3” correspond to CIB references 
models no. 1,2, and 3 shown in Fig. 1. 


HEGRA: preliminary 



Fig. 10. The intrinsic and absorbed 7 -ray spectra from 
1ES 1426+428 (see the text). 


6 The case of 1ES 1426+428 

Because of the relatively large redshift of 1ES 1426+428 
(z = 0.129), the TeV radiation from this extreme BL Lac 
object suffers severe IG-absorption. Therefore the discovery 
of TeV 7 -rays from this source by two independent (VER¬ 
ITAS and HEGRA) groups came, to a certain extent, as a 
pleasant surprise. The spectrum deformation factors calcu¬ 
lated for 3 different CIB models show that above 300 GeV 
the IG-absorption leads to significant steepening of the 7 -ray 
spectrum, but starting from 1 TeV to several TeV the spec¬ 
trum is deformed only slightly, although the suppression of 
the absolute flux may be as large as a factor of 100 (Fig. 9). 
This effect has a simple explanation discussed in Sec. 3. 

Remarkably, the preliminary spectrum of 1ES 1426+428 
derived from the HEGRA CT system data taken in 1999 and 
2000 (Gotting et al. , 2001), does show a steep slope at low 
energies with a tendency for a flattening above 1.5 TeV. The 
large statistical errors in flux measurements as well as the 
lack of information about the intrinsic 7 -ray spectrum pre¬ 
vent us from definite conclusions. Nevertheless, in Fig. 10 
we show the expected spectra of 7 -rays assuming a single 
power-law intrinsic spectrum with photon index identical to 
the index observed by BeppoSAX in the energy region up 
to 100 keV, T = 1.92 (Costamante et al., 2001). The solid 
curve represents the absorbed spectrum of 7 -rays calculated 
for the CIB model 1 and for the Hubble constant Hq = 
60 km/s/Mpc. In order to demonstrate the sensitivity of 
the intergalactic absorption effect to the ratio of the CIB flux 
and the Hubble constant, ucib/Hq (for the given redshift z, 
and the shape of the CIB spectrum, this ratio is the only pa¬ 
rameter which defines the optical depth), in Fig. 10 we show 
the 7 -ray spectra calculated for two other values of the ratio 
ucib/Hq which are 30% less (dashed curve) and 30% more 
(dot-dashed curve) compared to the the nominal value (solid 
curve). Hopefully, the future multi wavelength studies may 


allow us to distinguish between these three curves, thus we 
could be able to “measure” the absolute CIB flux at NIR, tak¬ 
ing into account that the current uncertainty of the the Hubble 
constant (presumably) does not exceed 30%. 

Independent of details, the pure fact of detection of TeV 
7 -rays from 1ES 1426+428 is an indicative of a need for re¬ 
vision of the current conceptual view of TeV blazars, accord¬ 
ing to which the synchrotron (X-ray) peak in the SED domi¬ 
nates over the inverse Compton (TeV) peak (see e.g. Fossati 
et al. 1998). Indeed, although the detected flux of 7 -rays is 
only few times 10 “ 12 erg/cm 2 s, corrected for IG-absorption 
this flux may well exceed 10 -10 erg/cm 2 s. For compari¬ 
son, the X-ray measurements (Costamante et al., 2001; T. 
Takahashi, private communication) show fluxes at the level 
of 10 -11 erg/cm 2 s. The TeV luminosity significantly ex¬ 
ceeds the X-ray luminosity, and apparently this should have 
an impact on the revision of the current TeV-blazar models. 

In particular, for the inverse Compton models (for review 
see Sikora and Madejski, 2001) this would require that the ra¬ 
diation of electrons occurs in the regime dominated by Comp¬ 
ton losses. Interestingly, this could serve as a good argu¬ 
ment for justification of the above assumption concerning 
the single power-law intrinsic TeV spectrum. Indeed, in this 
regime a pronounced feature (spectral flattening) formed in 
the steady-state electron distribution (caused by the struc¬ 
ture in the Klein-Nishina cross-section) is largely compen¬ 
sated by the same Klein-Nishina effect in the radiated 7 -ray 
spectrum; thus the E~ 2 type acceleration spectrum of elec¬ 
trons is (almost exactly) transfered to the 7 -ray spectrum (see 
e.g. Zdziarski and Krolik 1993). Of course, in this case we 
may expect, within the one-zone SSC model, quite hard syn¬ 
chrotron X-ray spectrum with spectral index a < 0.5. This 
actually does not contradict to the BeppoSAX data for which 
the broken power-law model gives very flat spectrum with an 






Fig. 11. Inverse Compton spectrum of cold unshocked ultrarela- 
tivistic jet with a Lorentz factor of bulk motion V = 3.33 x 10 7 . 
The radiation component associated with Comptonization of ambi¬ 
ent optical photons is shown by dashed curve I. The dashed curve II 
is the residual of the total TeV source emission (after subtraction of 
the unshocked wind component I), and formally could be attributed 
either to the Comptonization of the bulk motion on ambient FIR 
photons or to the radiation of blobs in shocked jet. The heavy solid 
line represents the superposition of these two components. Fits to 
the observed flux of Mkn 501 are shown by thin solid lines. Curve 1 
corresponds to the fit given by equation (b) in Fig. 3, and curve 2 - 
to the steepest possible spectrum above 17 TeV based on the recent 
reanalysis of Mkn 501 HEGRA data (Aharonian et al.. 2001a). The 
intrinsic (absorption-corrected) spectra of Mkn 501, corresponding 
to the fits to the observed points “1” and “2”, are shown by dotted 
lines “1'” and “2'”, respectively (from Aharonian et al., 2002). 


index above 10 keV, a-i «0.3-0.4 (Costamante et al., 2001). 

The hard intrinsic TeV spectra of 1ES 1426+428 with con¬ 
stant power-law slope and high L 7 /Lx > 10 ratio can be 
explained also by the proton-synchrotron model, assuming 
that a small, « 0.1 fraction of the proton-synchrotron TeV 7 - 
rays suffer internal absorption. If so, the hard X-ray emission 
could be attributed to synchrotron radiation of the secondary 
pair-produced electrons (Aharonian, 2000). The weak point 
of this model is in the requirement of very large, 100 G or so 
B-field and extremely high acceleration rate of protons. 

7 “IR background - TeV gamma-ray crisis” ? 

From the discussion of the previous section we may sum¬ 
marize that the basic temporal and spectral characteristics of 
TeV blazars can be described, at least qualitatively, by the the 
current leptonic, and perhaps also hadronic models. At the 
same time all these models would fail to explain the sharp 
pile-up which may appear at the end of the “reconstructed” 
spectrum of Mkn 501, if the reported FIR fluxes correctly de¬ 
scribe the level of CIB (see Fig. 3). Motivated by such a non¬ 
standard spectral shape, recently several extreme hypotheses 
have been proposed to overcome the “IR background - TeV 
gamma-ray crisis” (Protheroe and Meyer 2000). In particu¬ 


lar, Flarwit et al. (2000) suggested an interesting idea that 
the HEGRA highest energy events are due to Bose-Einstein 
condensations interacting with the air atmosphere. Another, 
even more dramatic hypothesis - violation of the Lorentz in¬ 
variance - has been proposed to solve this problem (e.g. Ki- 
fune, 1999; Stecker and Glashow, 2001). Two more propos¬ 
als could be added to the list of “exotic” solutions: (i) assum¬ 
ing that Mkn 501 is located at a distance <C 100 Mpc; or (ii) 
assuming that TeV 7 -rays from Mkn 501 are not direct repre¬ 
sentatives of primary radiation of the source, but are formed 
during the development of high energy electron-photon cas¬ 
cades in the intergalactic medium (Aharonian et al. 2001b). 
At first glance, both ideas seem relatively “neutral”. But in 
fact they do contain dramatic assumptions. While first hy¬ 
pothesis implies non-cosmological origin of the redshift of 
Mkn 501, the second hypothesis requires extremely low in¬ 
tergalactic magnetic field of order of 1 0 1 ,s G or less. 

Although very fascinating, the appeal for revisions of es¬ 
sentials of modern physics and astrophysics seems to be too 
premature. The nature of the FIR isotropic emission detected 
by COBE is not yet firmly established, and it is quite possible 
that the bulk of the reported flux, especially below 100 /iin, 
is a result of superposition of different local backgrounds. It 
should be noticed, however, that not only the flux at 60/zm, 
but also the CIB fluxes at longer, A > 100/rm wavelengths 
are responsible (albeit in a less distinct form) for the ap¬ 
pearance of the pile-up, unless we assume very specific CIB 
spectrum in the MIR-to-FIR transition region (Aharonian et 
al., 1999a; Coppi and Aharonian, 1999b; Guy et al., 2000; 
Renault et al., 2001). Therefore it is important to explore 
other possibilities for formation of pile-ups in the spectra of 
TeV blazars. Motivated by this, recently we have proposed 
a new, non-acceleration scenario which postulates that the 
> 10 TeV radiation of Mkn 501 is a result of the bulk mo¬ 
tion Comptonization of ambient low-frequency photons by a 
cold ultrarelativistic conical wind (Aharonian et al., 2002). 

The relativistically moving plasma outflows in forms of 
jets or winds, are common for many astrophysical phenom¬ 
ena on both galactic or extragalactic scales. Independent of 
the origin of these relativistic outflows, the concept of the 
jet seems to be the only successful approach to understand 
the complex features of nonthermal radiation of Blazars, Mi¬ 
croquasars and GRBs. The Lorentz factor of such outflows 
could be extremely large. In particular in the Crab Nebula 
the Lorentz-factor of the MHD wind is estimated between 
10 6 and 10' (Rees and Gun, 1974; Kennel and Coroniti, 
1984). Meszaros and Rees (1997) have shown that in the 
context of cosmological GRBs the magnetically dominated 
jet-like outflows from stellar mass black holes may attain 
extreme Lorentz factors exceeding 10 6 . The conventional 
Lorentz factors of jets in the inverse Compton models of 
7 -ray blazars, are rather modest, T ~ 10. However there 
are no apparent theoretical or observation arguments against 
the bulk motion with much larger Lorentz-factors (see e.g. 
discussion by Celotti et al., 1998). Because of existence of 
dense photons fields in the inner sub-parsec region of sources 
like Mkn 421 and Mkn 501, the Compton optical depth tq. 
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Fig. 12. Sketch of the two-stage 7 -ray production scenario in Mkn 501. At the first stage, the highest energy 7 -rays above 10 TeV with a 
sharp spectral form (“pile-up”) are produced due to Comptonization of the ambient optical photons by the cold ultrarelativistic jet with bulk 
motion Lorentz factor T ~ 10 7 — 10 8 . Low energy 7 -rays, E-, < 10 TeV can be also produced due to Comptonization of ambient FIR 
photons with broad spectral distribution. At the second stage, relatively low energy (E 1 < 10 TeV) 7 -rays are produced in the blobs moving 
with Lorentz factor of about L ~ 10 in accordance with the standard SSC model (from Aharonian et al., 2002). 


could be as large as 1. Obviously, for our model the most fa¬ 
vorable value for rc lies between 0.1 (in order to avoid huge 
energy requirements to the outflow) and 1 (in order to avoid 
the Compton drag). In particular, we have to assume that in 
the proximity of the black hole the outflow should be Poynt- 
ing flux dominated, and only at large distances from the cen¬ 
tral object, where the photon density is significantly reduced, 
the major part of the electromagnetic energy is transfered to 
the kinetic energy of bulk motion. 

Due to the extremely large Lorentz factors exceeding 10 7 , 
the Compton scattering on the ambient NIR/optical photons 
with energy more than 1 eV proceeds in deep Klein-Nishina 
regime; therefore the 7 -radiation should have a very narrow 
distribution with energy E ss E e = m e c 2 Y. Meanwhile, 
the IC scattering on ambient far IR photons surrounding the 
central source, still takes place in the Thomson regime, and 
thus results in a smooth broad-band spectrum. 

Fig. 11 demonstrates that the overall absorption-corrected 
(for the CIB at FIR close to the extreme CIB model no.3) 
spectrum of Mkn 501 can be satisfactorily explained in the 
terms of the inverse Compton emission of the cold jet with 
bulk Lorentz factor F = 3.33 x 10 7 , assuming an ambi¬ 
ent radiation field with a narrow (Planckian) type radiation 
with temperature kT = 2 eV (dashed curve I). The dashed 
curve II which formally is the residual from the subtraction 
of the unshocked wind component radiation from the intrin¬ 
sic (reconstructed) TeV emission (solid heavy line), can be 
attributed to the bulk Comptonization on ambient far IR pho¬ 
tons produced e.g. by cold clouds surrounding the source. 
Alternatively, the “residual component” (curve II) can be re¬ 
ferred to the SSC (or any other) radiation component of blobs 
in shocked jet. This two-stage ( pre-shock plus post shock ) 
scenario of formation of TeV 7 -ray emission is schematically 
illustrated in Fig. 12 

The possibility to disentangle the multi-TeV emission with 
a characteristic sharp pile-up at the very end of the spectrum, 
Ej ~ m e c 2 r, from the sub-10 TeV emission associated with 
the shocked structures (e.g. blobs) in the jet, not only solves 
the possible “IR background - TeV gamma-ray crisis” but 
also allows more relaxed parameter space for interpretation 


of X-rays and the remaining “low” energy (< 10 TeV) 7 - 
rays within the conventional SSC scenario. Consequently, 
this offers more options for interpretation of X-ray/gamma- 
ray correlations. If the overall TeV radiation of Mkn 501 
indeed consists of two, unshocked and shocked jet emission 
components, we may expect essentially different time behav¬ 
iors of these components. In particular, the “unshocked jet” 
(> 10 TeV) component should arrive earlier than the SSC 
components consisting of synchrotron X- and sub-10 TeV 
IC 7 -rays. Therefore, an important test of the suggested 
two stage scenario of the TeV radiation of jets would be 
the search for correlations (or lack of such correlations) of 
> 10 TeV radiation with both the low energy (e.g. 1-3 
TeV) 7 -rays and synchrotron X-rays. The low statistics of 
(heavily absorbed) 7 -rays above 10 TeV makes the search 
for such correlations rather difficult, and requires ground- 
based instruments with huge, ^> 0.1 km 2 detection areas in 
this energy domain. The new generation IACT arrays like 
CANGAROO-3, H.E.S.S. and VERITAS should be able to 
perform such correlation studies. 

It should be noticed, however, that the study of the sig¬ 
natures of IG-absorption of 7 -rays from Mkn 421 and Mkn 
501 does not contain adequate information about CIB at A > 
IO// 111 . This information is essentially lost due to the heav¬ 
ily absorbed 7 -rays above 20 TeV. This objective is rather 
contingent on discovery of nearby extragalactic TeV sources 
within « 50 Mpc, dictated by the condition that the mean 
free path of 7 -rays at these energies should not significantly 
exceed the distance to the source (see Fig. 2). Despite the 
lack of very close blazars, other potential extragalactic 7 -ray 
sources like the nearby radiogalaxies Centaurus A and M 87, 
and perhaps also the startburst galaxies Arp 220, M 82 and 
NGC 253 may (hopefully) provide us with multi-TeV 7 -rays 
for such important studies. 

Yet, the distant TeV sources may offer an additional, and 
perhaps even more informative channel for derivation of CIB 
fluxes at FIR wavelengths. Such a possibility is based on the 
study of angular and spectral characteristics of hypothetical 
giant Pair Halos expected around powerful nonthermal ex¬ 
tragalactic sources (Aharonian et al., 1994). 
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The figure is available as a separate file 

Fig. 13. Illustration of the process of formation of Pair Halos. 


8 Pair Halos and CIB 

Because of IG-absorption, the bulk of primary TeV 7-rays 
from distant extragalactic sources is not visible. The energy 
of these photons however is not lost - the electromagnetic 
cascade initiated by the absorbed 7-rays can be observed ! 
The cascade in the IG medium (e.g. Berezinsky and Kudri¬ 
avtsev, 1990; Protheroe and Stanev, 1993; Coppi and Aha- 
ronian, 1997) is supported by 2 processes: (i) by inverse 
Compton scattering on the photons of 2.7 K MBR and (ii) by 
photon-photon pair production on CIB (for energies above 
100 TeV - also on 2.7 K MBR). When the magnetic field 
near the source is sufficiently large, B > 10~ 9 G (generally 
on the scales of tens of Mpc, but if the spectrum of primary 7- 
rays extends into PeV region - within tens of kpc), i.e. when 
the mean free path of cascade electrons (at least at the ini¬ 
tial stages) is large compared to Larmor radii, the cascade 
will generate, even for a highly beamed primary source, an 
extended isotropic Pair Halo (see Fig. 13). While the forma¬ 
tion of Pair Halos is unavoidable for all extragalactic objects 
with spectra extending beyond 10 GeV, only relatively com¬ 
pact Pair Halos formed around (multi) TeV sources could 
be detectable. The superposition of the cascade radiation of 
all Pair Halos formed in the Universe may contribute sig¬ 
nificantly to the diffuse GeV 7-ray background (Coppi and 
Aharonian, 1997). 

The halo radiation can be distinguished by its characteris¬ 
tic variation in spectrum and intensity with angular distance 
from the central source, which weakly depends on the details 
of the source model, in particular on the energy spectrum and 
the geometry of the primary source, especially when a hard 
primary 7-ray spectrum extends to PeV energies. 

The energy E 1 and the arrival angle 8 of the detected pho¬ 
ton from a source at a distance d are determined, on average, 
by the energy of the parent electron E e « (75 7 /4/cTo) 1 / 2 m e c 2 
(To = 2.7 K), and the distance l from the central source 
where this electron is produced. Because the mean free path 
of 7-rays A decreases with energy (see Fig. 2), the energy 
and site of production of the secondary electron are defined, 
on average, by the “last 7-CIB interaction”, E e « E lt o/2, 
therefore the energy E-y and arrival angle 8 of the detected 7- 
rays contain information about their “grandparents” - 7-rays 


with energy E 7i o ~ 34(T 7 /1 TeV) 1 / 2 TeV, and mean free 
path A(T 7 > o) = d ■ 8 . Correspondingly, the measured mean 
free path A(T 7 j o) should tell us about the density of CIB at 
A ~ 50(£7 7 /l TeV) 1 / 2 /an. Thus, detections of Pair Halos 
surrounding Mkn 421 and Mkn 501 at 1 TeV should provide 
information on CIB around 50 pm, provided that the primary 
spectrum extends beyond 30 TeV. 

Recently the HEGRA collaboration made the first attempt 
to probe the existence of a Pair Halo in direction of Mkn 501, 
using the data base accumulated during the period of obser¬ 
vations from 1997 to 1999 (Aharonian et ah, 2001c). The de¬ 
rived limit on the halo flux Jph(> ITeV) ~ 10 _ 8 ph/cm 2 s sr 
within 1 ° is unfortunately above the theoretical expectations 
based on the energy budget constraints. Indeed, Mkn 501 
actually is not an ideal candidate for a halo source given its 
relatively modest luminosity and the fact that the preferable 
distance to observe halo is in the range of several 100 Mpc to 

1 Gpc; consequently, the most effective searches can be per¬ 
formed by detectors with energy threshold around 100 GeV 
(Aharonian et ah, 1995). In this regard, the BL Lac objects 
1ES 1426+428 at 2 = 0.129 and PKS 2155-30 at z = 0.116 
with 7 -ray luminosities exceeding the “quiescent” TeV lumi¬ 
nosity of Mkn 501 al least by a factor of 10, seem to be more 
suitable objects for producing detectable Pair Halos. The 
forthcoming H.E.S.S. array of 10m IACTs equipped with op¬ 
timal for this purpose 5° FoV high resolution cameras, will 
be the first instrument for a sensitive search for Pair Halos at 

2 ~ 0.1-0.2. 

Generally, the blazars are are relatively weak sources - 
they seem very bright because of the Doppler boosting (pro¬ 
portional to ()j l > 10 3 ). The distant powerful radiogalaxies 
and quasars with large-scale jets (as potential accelerators of 
ultrahigh energy cosmic rays) may appear more promising 
sites for production of detectable - strong and compact - Pair 
Halos (Aharonian, 2002). The detection and identification of 
such halos from cosmologically distant objects beyond z = 1 
is possible at energies of 7 -rays of about 10 GeV. While the 
spectral cutoffs in the energy spectra of the direct and halo 
emission components are caused by absorption in the diffuse 
UV background radiation, the angular size of the halo con¬ 
tains information about the CIB at several to 10 pm wave¬ 
lengths at the epoch corresponding to the redshift of the cen¬ 
tral source 2. This is a unique (not available by other means) 
channel which provides us with model-independent informa¬ 
tion about the cosmological evolution of CIB. Such studies 
can be effectively performed by GLAST (e.g. Gehrels and 
Michelson 1999) and by future sub-10 GeV threshold high- 
altitude IACT arrays like 5 @5 (Aharonian et al. 2001). 

9 Summary 

The energy-dependent mean free path of 7 -rays in the in- 
tergalactic medium at TeV energies does not exceed several 
100 Mpc. Therefore VHE 7 -rays from TeV blazars arrive 
with significantly distorted spectra. Our limited knowledge 
of CIB results in large uncertainties in the reconstructed (cor- 
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rected for intergalactic absorption) intrinsic 7 -ray spectra. 
Consequently, in spite of good quality of spectral measure¬ 
ments of two strongest TeV blazars, Mkn 421 and Mkn 501, 
nowadays we are faced with a challenge - our understanding 
of radiation processes in relativistic jets is neither complete 
nor conclusive. Nevertheless, the gamma-ray astronomers 
believe that eventually (soon ?) they will learn to identify 
confidently the radiation mechanisms, to fix/constrain the rel¬ 
evant model parameter space, and calculate robustly the in¬ 
trinsic 7 -ray spectra based on the multiwavelength studies of 
spectral and temporal characteristics of blazars obtained si¬ 
multaneously at X-ray and TeV bands on sub-hour timescales. 
Then it will be possible to estimate unambiguously the effect 
of the intergalactic 7 -ray absorption, and thus to infer robust 
information about the CIB fluxes. Moreover, the studies of 
angular and spectral properties of giant Pair Halos formed 
around powerful extragalactic multi-TeV sources may pro¬ 
vide us with a unique tool for derivation of the spectra and 
absolute fluxes of CIB at different cosmological epochs z, 
and thus to probe the evolution of galaxies in past. 
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